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Abstract 

The subject of this work is the patrolling of an environment with the aid of a team of autonomous agents. 
We consider both the design of open-loop trajectories with optimal properties, and of distributed control laws 
—i converging to optimal trajectories. As performance criteria, the refresh time and the latency are considered, i.e., 

respectively, time gap between any two visits of the same region, and the time necessary to inform every agent 
about an event occurred in the environment. We associate a graph with the environment, and we study separately 
the case of a chain, tree, and cyclic graph. For the case of chain graph, we first describe a minimum refresh time 
and latency team trajectory, and we propose a polynomial time algorithm for its computation. Then, we describe a 
distributed procedure that steers the robots toward an optimal trajectory. For the case of tree graph, a polynomial 
time algorithm is developed for the minimum refresh time problem, under the technical assumption of a constant 
number of robots involved in the patrolling task. Finally, we show that the design of a minimum refresh time 
trajectory for a cyclic graph is NP-hard, and we develop a constant factor approximation algorithm. 
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I. Introduction 

The recent development in the autonomy and the capabilities of mobile robots greatly increases the 
number of applications suitable for a team of autonomous agents. Particular interest has been received 
by those tasks requiring continual execution, such as the monitoring of oil spills [lj, the detection of 
forest fires 0, the track of border changes 0, and the patrol (surveillance) of an environment 0. 
The surveillance of an area of interest requires the robots to continuously and repeatedly travel the 
environment, and the challenging problem consists in scheduling the robots trajectories so as to optimize 
^ a certain performance criteria. The reader familiar with network location, multiple traveling salesman, or 
graph exploration problems may observe a close connection with the patrolling problem we address, e.g., 

, \ see 0, 0, 0. It is worth noting, however, that these classical optimization problems do not capture the 

O repetitive, and hence dynamic, aspect of the patrolling problem, nor the synchronization issues that arise 
when a timing among the visits of certain zones is required. 

A precise formulation of the patrolling problem requires the characterization of the robots capabilities, 
of the environment to be patrolled, and of the performance criteria. In this work, we assume the robots 



^ to be identical and capable of sensing and communicating within a certain spatial range, and of moving 
jf3 according to a first order integrator dynamics with bounded speed. We represent the environment as a 
graph, in which the vertices correspond to physical and strategically important locations, and in which the 
edges denote the possibility of moving and communicating between locations. We assume that, when a 
robot is placed at each of the graph vertices, the union of the sensor footprints provides complete sensor 
coverage of the environment. Regarding the performance criteria of a patrolling trajectory, we consider (i) 
the time gap between any two visits of the same region, called refresh time, and (ii) the time needed to 
inform the team of robots about an event occurred in the environment, called latency. Loosely speaking, 
refresh time and latency reflect the effectiveness of a patrolling team in detecting events in the environment 
and in organizing remedial actions. For both the refresh time and latency optimization problem, we focus 
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on the worst case analysis, even though the average refresh time and the average latency cases are also 
of interest. Notice that for the latency to be finite, the motion of the robots needs to be synchronized. For 
instance, if two robots are allowed to communicate only when they simultaneously occupy two adjacent 
vertices of the graph, then they need to visit those vertices at the same time in a finite latency trajectory. 

The patrolling problem is receiving increasing attention because of its fundamental importance in many 
security applications, e.g., see (H, fl9), ifTOl . ifPTfl. [[T2]|. Although many solutions have been proposed, the 
problem of designing minimum refresh time and latency team trajectories for a general environment is, 
to date, an open problem. Almost all traditional approaches rely on space decomposition, and traveling 
salesperson tour computation [fT3l . In lfl4ll an empirical evaluation of existing patrolling heuristics is 
performed. In [TT5l two classes of strategies are presented, namely the cyclic- and the partition-based 
strategy. In the cyclic-based strategy, the robots compute a closed route through the viewpoints, and travel 
repeatedly such route at maximum speed. Clearly, in the case of a single robot, if the tour is the shortest 
possible, then the cyclic-based strategy performs optimally with respect to the refresh time and latency 
criteria. In the partition-based strategy, the viewpoints are partitioned into m subsets, being m cardinality 
of the team, and each robot is responsible for a different set of viewpoints. To be more precise, each robot 
computes a closed tour visiting the viewpoints it is responsible for, and it repeatedly moves along such 
tour at maximum speed. Still in [fT5ll . the two classes of strategies are compared, and it is qualitatively 
shown that cyclic-based strategies are to be preferred whenever the ratio of the longest to the shortest 
edge of the graph describing the environment is small, while, otherwise, partition-based policies exhibit 
better performance. In fl4[ and Q, an efficient and distributed solution to the perimeter patrolling problem 
for robots with zero communication range is proposed. By means of some graph partitioning and graph 
routing techniques, we extend the results along these directions, e.g., by considering the case of a nonzero 
communication range for the perimeter patrolling, and by characterizing optimal strategies for different 
environment topologies. An important variant of the patrolling problem is known as persistent surveillance, 
e.g., see IfToll . Differently to our setup, a dynamically changing environment is considered for the persistent 
surveillance problem, and performance guarantees are offered only under a certain assumption on the rate 
of change of the regions to be visited. 

It is worth mentioning that a different approach to the design of patrolling trajectories relies on the use 
of pebbles or chemical traces to mark visited regions, e.g., see [TTTll . lfT8l . [fT9l . EOl . These techniques, 
although effective even without a global representation of the environment, and with severe communication 
constraints, do not explicitly deal with the optimality of the patrolling trajectories, and they represent 
therefore a complementary area of research with respect to this work. 

The main contributions of this work are as follows. We introduce and mathematically formalize the 
concept of refresh time and latency of a team trajectory, and we formally state the patrolling optimization 
problem. We propose a procedure to build a graph (roadmap) to represent the topological structure of the 
area to be patrolled, and we study separately the case of a chain, tree, and cyclic (not acyclic) graph. We 
exhaustively discuss the case of a chain roadmap. First, we characterize a family of minimum refresh time 
and latency team trajectories, which can be computed by optimally partitioning the chain graph among 
the robots. Second, we derive a centralized polynomial time algorithm to compute an optimal partition, 
and, ultimately, to design an optimal team trajectory. Our partitioning procedure is based upon a bisection 
method, and it is also amenable to distributed implementation. Third, we develop a distributed procedure 
for the robots to converge and synchronize along an optimal trajectory, so as to minimize the refresh time 
and latency criteria. Fourth and finally, we test the robustness of our methods through a simulation study. 
When the roadmap has a tree or cyclic structure, we focus on the refresh time optimization problem, and 
we do not consider the latency optimization nor the design of distributed algorithms. For the case of a tree 
roadmap, we reduce the minimum refresh time patrolling problem to a known graph optimization problem. 
We show that the computational complexity of the minimum refresh time patrolling problem is polynomial 
in the number of vertices of the roadmap, and, under the assumption of a fixed and finite number of robots, 
we identify a polynomial time centralized algorithm to compute a minimum refresh time team trajectory. 
For the case of a cyclic roadmap, we show that the patrolling problem is an NP-hard optimization problem. 
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We propose two approximation algorithms, and we characterize their performance. The first approximate 
solution is extremely easy to compute, but its performance depends upon the ratio between the longest 
and the shortest edge in the graph representing the environment. The second approximation algorithm is 
based on a polynomial time path-covering procedure, and it allows us to compute a team trajectory whose 



refresh time is within a factor of 8 from the minimum refresh time for the given environment (cf. Fig. 13 
for an example). To the best of our knowledge, this algorithm is the first constant factor approximation 
algorithm for the NP-hard minimum refresh time patrolling problem. 

A preliminary version of this work appeared in [f2TTl . With respect to the latter manuscripts, in this 
current work we introduce and solve the latency optimization problem, we perform a numerical study to 
analyze the robustness of our algorithmic procedures, and we improve the presentation of the results on 
the refresh time optimization problem. 

The rest of the paper is organized as follows. The notation and the problem under consideration are in 
Section [EEJ where we also show that the patrolling problem is, generally, computationally hard. Section 



III IV] and [V] contain our results for the patrolling of a chain environment. We characterize a minimum 



refresh time and latency team trajectory, and we derive a centralized and a decentralized algorithm for its 



computation. In Section VI we perform a simulation study to show some robustness and reconfigurability 



properties of our distributed procedure. Section VII-A contains our results for the patrolling of a tree 



environment. We describe a minimum refresh time team trajectory on a tree roadmap, and we characterize 



the complexity of computing an optimal solution. Section VII-B deals with the general case of cyclic 
environment, and it contains our approximation procedures. Our conclusion and final remarks are in 
Section [VTTT1 

II. Robotic model and preliminary concepts 
A. Robots on roadmaps with sensor coverage and communication connectivity 

We will be using the standard motion planning notation, and we refer the reader to Il22l for a com- 
prehensive treatment of the subject. We are given a team of m > identical robots, capable of sensing, 
communicating, and moving in a connected environment. We make the following combined assumptions 
on the robot capabilities and on the environment to be patrolled. 

Regarding sensing, we assume that the environment can be completely covered by simultaneously 
placing a robot at each of a set of n > m viewpoints in the configuration space. In other words, if m = n 
robots were available and placed at the n viewpoints, then the union of the sensors footprint of the robots 
would provide complete sensor coverage of the environment. We assume that each viewpoint is required 
for complete sensor coverage. Finally, we assume n > m so that at least one robot needs to visit more 
viewpoints for the entire environment to be monitored over time. 

Regarding communication, we associate an undirected graph G with the environment, whose vertices 
are the n viewpoints, and in which there is an edge between two vertices if two robots placed at those 
viewpoints are able to communicate to each other. We assume that G is connected. In what follows 
we design cooperative patrolling algorithms with sporadic communication, in the sense that two robots 
are required to communicate only when they occupy adjacent vertices. The occurrence of additional 
communication links can be easily incorporated into the algorithms and cannot decrease their performance. 

Regarding motion, we assume that the robots are holonomic, i.e., they are modeled as first order 
integrators and move at most at unit speed. Additionally, we turn the graph G into a robotic roadmap [|22| 
and a metric weighted graph as follows: to each pair of viewpoints that are neighbors in G, we associate 
a unique path connecting them. We select these paths so that the set of path lengths, adopted as edge 
weights, verify the triangle inequality. (For example, the shortest paths between viewpoints constitute a 
suitable choice.) We assume that each robot remains always on the roadmap. 

In summary, the combined assumptions on the robot capabilities and on the environment are that: the 
vertices of G provide complete sensor coverage of the environment and each edge of G corresponds to 
both a communication edge and a motion path. Hence, we refer to G as a roadmap with sensor coverage 
and communication connectivity. 



Fig. 1. A polygonal environment and its associated roadmap. The viewpoints are chosen in a way that the environment is completely 
covered by placing a robot at each of the 12 viewpoints. The edges of the roadmap denote the possibility for the 3 robots of both moving 
and communicating between pair of connected viewpoints. The weight of each edge corresponds to its length. 



B. On suitable roadmaps and robots capabilities 

The problem of constructing a roadmap from an environment is here discussed. 

Example 1 (Roadmap computation with omnidirectional sensing and communication in known environ- 
ments). Assume that the robots are holonomic vehicles moving at bounded speed, and equipped with an 
omnidirectional sensing device, and a line-of-sight communication device. If a map of the environment is 
available, then a valid roadmap is obtained by solving an Art Gallery Problem with Connectivity / [23l/ . 
A solution to the Art Gallery Problem with Connectivity is a set of locations, called guards, with the 
following two properties: each point in the environment is visible by at least one guard and the visibility 
graph of the guards is connected. An example roadmap is given in Fig. [7] (A distributed sensor-based 
algorithm for the Art Gallery Problem with Connectivity is given in K2M .) 



Example 2 (Roadmap computation for general robots in an unknown environment). If the environment is 
unknown, then an exploration algorithm should be used to obtain a representation of the environment, see 
for example [18], / l25l/ . / l26l/ . While or after exploring the environment, a robot should select viewpoints 
to provide full sensor coverage. By construction, the existence of a path between two viewpoints is 
automatically guaranteed by the exploration algorithm. Moreover, if communication between two path- 
connected viewpoints (vi, Vj) is not guaranteed, then the graph may be augmented with additional 
viewpoints along the path connecting and vj. The roadmap resulting from these steps features sensor 
coverage and communication connectivity. 

The following remarks are in order. First, a roadmap representing the environment is in general not 
unique, and the patrolling performance depends upon the particular choice. In this work, we do not 
address the problem of choosing the roadmap that yields optimal performance. Instead, we present efficient 
algorithms for the patrolling problem, which can also be used to compare different roadmaps on the basis 
of the corresponding patrolling performance. Second, for the implementation of our algorithms, a robot 
does not need to travel exactly along the roadmap. Indeed, a robot only needs to arrive sufficiently close 
to a viewpoint, or to be able to circle around it to provide sensor coverage. A related example is in 



Section VI where the arrival delay of a robot can be interpreted as the uncertainty in the motion of the 
robots. Third, the global knowledge of the roadmap may not be required for the execution of a patrolling 
trajectory. Indeed, in general, each robot only visits a subpart of the roadmap. Fourth and finally, collisions 
are prevented by letting two robots exchange their roles every time they are about to collide. Indeed, since 
no robot is assigned to a specific viewpoint, if robot i is about to collide with robot j at time T, then, 
without affecting the performance of the team trajectory, collision is avoided by redefining the 2-th and 
j-th trajectory as Xi(t) = Xj(t) and Xj(t) = Xi(t) for t > T; see the notion of order invariant trajectory 
below. Communication or sensing devices can be used to detect and prevent possible collisions. 
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C. Refresh time of team trajectories: definitions and a preliminary complexity result 

For a team of m robots with specific capabilities, let G = (V, E) be a robotic roadmap with sensor 
coverage and communication connectivity. A team trajectory X is an array of m continuous and piecewise- 
differentiable trajectories xi(t),... ,x m (t) defined by the motion of the robots on the roadmap G, i.e., 
Xi : [0, Tf] I—)- G, for i G {1,. . . ,m}, where T f G R is a time horizon of interest, much larger than the 
time required by a single robot to visit all viewpoints. We say that a viewpoint v G V is visited at time t 
by robot i if Xi(t) = v. We define the refresh time of a team trajectory X, in short RT(X), as the longest 
time interval between any two consecutive visits of any viewpoint, i.e., 

RT(X) = max max t 2 — ti 

vev (t 1 ,i 2 )en(t),x) 

where Q(v,X) = {(ti,t 2 ) G [0,T { ] 2 ,ti < t 2 , | x^t) ^ v, Vi G {1, . . . ,m},h < t < t 2 }. 

Remark 1 (Existence of a minimum). We claim that there exists a team trajectory with minimum refresh 
time and prove it as follows. Without loss of generality, we restrict our attention to team trajectories in 
which each robot moves at maximum speed along the edges and stops for certain durations, possibly 
equal to zero, at the viewpoints. Thus, a team trajectory can be described as a tuple of the form (S, A) = 
{(Si, Ai), . . . , (S m , A m )}, where Si contains the sequence of viewpoints visited by robot i, and Aj contains 
the waiting times at the visited vertices. Notice that the time horizon Tf is finite, the length of each edge is 
positive, the number of vertices is finite, and the speed of the robots is bounded. It follows that the length 
of each sequence Si is finite, and, therefore, each Si takes value in a finite set. Now, for each possible 
sequence of visited vertices, the refresh time is a continuous function of only the waiting times, and each 
waiting time lies in the compact interval [0, Tf]. Because any continuous function defined over a compact 
region admits a point of minimum value, we conclude that there exists a team trajectory with minimum 
refresh time. 

Problem 1 (Team refresh time). Given a roadmap and a team of robots, find a minimum refresh time 
team trajectory. 



In Section IV we present a different optimization problem, which deals with the possibility for a robot 
to communicate, possibly with multiple hops, with every other robot in the team. We now conclude this 
section with our first result on the computational complexity of the Team refresh time problem. For a 
detailed discussion of the main aspects of the computational complexity theory, we refer the interested 
reader to fl27). Recall that an optimization problem is said to be NP-hard if it is, informally, as hard as 
the hardest problem in the class NP, for which no polynomial time algorithm is known to compute an 
optimal solution. 

Theorem II.l (Computational complexity). The Team refresh time problem is NP-hard. 

Proof: This statement can be shown by reduction from the Traveling Salesman problem Il27il . In fact, 
if m = 1, since the speed of the robots is bounded, then a minimum refresh time trajectory consists of 
moving the robot at maximum speed along a shortest closed tour visiting the viewpoints. The problem of 
finding a shortest tour through a set of points in the plane, also known as Traveling salesman problem, is 
an NP-hard problem ll27Tl . Hence, by restriction, the Team refresh time problem is also NP-hard. ■ 



Following Theorem II. 1[ the minimum refresh time optimization problem is generally computationally 



hard. In this work, we first identify two roadmap structures for which there exists an efficient solution 
to the Team refresh time problem, and then we describe two approximation algorithms to deal with the 
general case. 

III. Minimum refresh time team trajectory on a chain roadmap 

We characterize in this section an optimal refresh time team trajectory when the roadmap associated 
with the environment has a chain structure. 
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A. Open loop team trajectory characterization 

Let Ni denote the neighbor set of the vertex i, and let |A^| denote the degree of i, i.e., the cardinality 
of the set N^. A chain roadmap is an undirected, connected, and acyclic roadmap, in which every vertex 
has degree two, except for two vertices which have degree one. Without losing generality, we assume 
that the n vertices are ordered in a way that N\ = {2}, N n = {n — 1}, and Ni = {i — l,i + 1} for 
each i E {2, . . . ,n — 1}. We define a relative order of the robots according to their position on the 
roadmap. A team trajectory is order invariant if the order of the robots does not change with time, i.e., 
if Xi(t) < Xi + i(t) for each i E {1, . . . ,m — 1} and for every instant t E [0, Tf], where Xi(t) denotes the 
distance at time t on the roadmap from the first vertex of the chain to the position of the i-ih. robot. 

Proposition III.l (Order invariant team trajectory). Let X be a team trajectory. There exists an order 
invariant team trajectory X such that RT(X) = RT(X). 

Proof: Let X be a team trajectory, and consider the permutation matrix P(t), that keeps track of the 
order of the robots at time t, i.e., such that the (z,j)-th entry of P(t) is 1 if, at time t, the 2-th robot 
occupies the j-th position in the chain of robots, and it is otherwise. Since X is continuous, anytime the 
function P(t) is discontinuous, the positions of the robots directly involved in the permutation overlap. 
Therefore, the order invariant team trajectory X = P~ 1 (t)X(t) is a feasible team trajectory, and it holds 
RT(X) = RT(X). " ■ 

Let Vi C V be the set of viewpoints visited over time by the agent i with the trajectory Xi, and let the 
image of the team trajectory X be the set {Vi, . . . , V^}. Notice that different team trajectories may have 
the same image. Throughout the paper, let Zj = min„ e y. v, r\j = max ue y v, and di — r» — k. Finally, let 
RT* = minx RT(X). A team trajectory is non-overlapping if V; fl V} = for all i ^ j. 

Proposition IIL2 (Non-overlapping team trajectory). Given a chain roadmap, there exists an order 
invariant and non- overlapping team trajectory with refresh time RT*. 

Proof: Let X* be a minimum refresh time team trajectory, and let X be the order invariant team 



trajectory obtained from X* as in Proposition III.l Clearly RT(X) = RT*. Let {Vi, . . . , V m } be the image 
of X, and note that V = U^Vj. Consider the partition of V defined as 

Vy = Vi, 

V i = V i \^y i , ze{2,...,m}. 

For every nonempty Vi, let ^ = min„ 6 y. v, f» = max^y. v, and Jj = fj — Ij. Note that, by construction, 
the viewpoint ^ is visited by the robot i and, possibly, by the robots j > i. Also, because X is order 
invariant, we have X{(t) < Xj(t). It follows that RT(X) > 2maxjJj. Consider now the team trajectory 
X with image {Vi,...,^}, and assume that the robots sweep periodically at maximum speed their 
segment. Then RT(X) = 2maxjrfj, so that X is an order invariant and non-overlapping team trajectory 
with minimum refresh time. ■ 
Given a chain graph on the viewpoints V, let IT m = {tti, . . . , Tc m } be an m-partition of V, i.e., 7Ti, . . . , 7r m 
is a collection of subsets of V such that TiiHTij = whenever i ^ j, and V = UI^=i 71 i- Additionally, let the 
dimension of the partition Ii m equal the longest distance between any two viewpoints in the same cluster, 
i.e., dim(n m ) = max,g {1| ,„ im } (max„ g7ri v - min„ e7r . v), where max„ 6iri t; - min ve7Ti v = if ^ = 0. 



Following Proposition III.2, there exists a minimum refresh time team trajectory whose image coincide 
with an m-partition of V. We now show that the minimum refresh time equals twice the dimension of an 
optimal m-partition. 

Theorem III.3 (Minimum refresh time). Let G be a chain roadmap, and let m be the number of robots. 
Then RT* = 2 minn m dim(II m ). 



Proof: As a consequence of Propositions |III.l and III.2 there exists a minimum refresh time team 



trajectory whose image coincides with an m-partition Il m . Since each robot is assigned a different cluster, 
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Trajectory 1: Minimum refresh time trajectory on a chain roadmap (i-th robot) 

Input : U := min l)6 y- i v, Vi := max„ 6 Vi v, di := Ti — U\ 
Require : an optimal partition of the chain graph; 

i: Xi(t) := h for t := 0, 2d>, 4dj, . . . ; 
2: Xi(t) := Ti for t :— di, 3di, 5d iy . . . ; 
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(a) 



(b) 

Fig. 2. A 4-partition with minimum (maximum) dimension is in Fig. |2(a)| The dimension of this partition is 117 — v§. A 4-partition with 
minimum (average) dimension is obtained by removing 3 longest edges and it is reported in Fig. |2(b)| The dimension of this partition is 

V10 — V7 > V7 — V(,. 



and the speed of the robots is bounded by 1, we have RT* > 2dim(il m ). Consider a team trajectory 
X in which each robot continuously sweeps at maximum speed the cluster it is assigned to. Clearly, 
RT(X) = RT* = 2dim(n m ). " ■ 

We have shown that a minimum refresh time trajectory consists of letting the robot sweep at maximum 
speed a part of the chain graph. Such a trajectory is more formally described for the i-th robot in Trajectory 
[T] where we only characterize the instants of time at which robot i changes its velocity vector, and we 
assume that it moves at maximum speed otherwise. 

Remark 2 (Average partition). By removing the longest edges in the chain the average length of the 
clusters is minimized. In general, such partition does not minimize the dimension of the m-partition, and 
hence it is not optimal in our sense. An example is in Fig. [2] 

B. Optimal m-partition centralized computation 

In the remaining part of the section we describe an algorithm to compute an optimal m-partition. For 
a set of viewpoints V, we call left-induced partition of length p E IR>o the partition U p = {ni} defined 
recursively as (cf. Fig. |3(a)| ) 

TTi = {V e V I CI, < V < di + p}, i = 1, . . . , (1) 



where 



min{t> G V \ v > cij_i + p}, 



The cardinality |il p | corresponds to the integer j such that {v 6 V \ v > aj + p] = 0. Observe that the 



function p i-)- |I1 P | is monotone, non-increasing, and right-continuous (cf. Fig. 3(b)). Let . . . ,p n _i} G 



i>n— 1 
">0 



be the discontinuity points of the function p i-)- |IF|, then, for k E {1, . . . ,n — 1}, 



(2) 



\Tl p \<k, if p>p k , 
|rp| > k, if p < p k . 

Note that two or more discontinuity points of |n p | may coincide, so that the function |n p | may not assume 



all the values in the set {!,... ,n}, e.g., in Fig. 3(b) the value |il p | is never equal to 9 
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Fig. 3. In 3(a) the left-induced partition of length p4 and p$, with ps < pi, for the chain roadmap with vertices {vi, . . . , vio}. The 
cardinalities are | II P4 1 = 4 and |II P5 | = 5, respectively. In |3(b)| the cardinality |n p | is plotted as a function of the length p. Notice that, 
because v 2 — vi — Vio — vg, the function |n p j does not assume the value 9. 



Algorithm 2: Optimal left-induced m-partition 



Input 

Set 



Viewpoints Vi , . 



. , v„, Number of robots m < n, Tolerance < e < — ; 

n ._ (<t+6) . 

f ■ 2 ' 



while (b — a) > 2e do 

n p := left-induced({ni, 
if |n p | > m then 

| a:=p,p:=^; 
else 

[ IT := n p , b:= p, p := 



6-. return IT* 



a + b . 
2 ' 



Theorem III.4 (Optimal m-partition). Let G be a chain roadmap. Let Ii Ti 
let Il p be the left-induced partition of length p of G. Then 



be an m-partition of G, and 



min dim(Il m ) = min{p G 



^>o 



|n p | < m}. 



Proof: Let U m be an m-partition, and let U p = {n^, . . . , 7r£} be the left induced partition of length 
p of a chain roadmap G. Let p* = min Hm dim(Il m ). We want to show that p* is one of the discontinuity 
points of the function |IP|, i.e., that p* verifies the conditions (|2]). 

By contradiction, if p < p* and |I1 P | < m, then an m-partition with dimension smaller than the optimal 
would exists. Therefore we have |n p | > m if p < p* . 

Suppose now that p > p* , and let = {nl, . . . ,7rj^} be an m-partition with minimum dimension. 
Notice that |7rf| > |7r*|, because the cluster 7rf contains all the viewpoints within distance p from v±, 
and hence also within distance p*. It follows that max 71^ > max7r*, and also that minvr^ 1 > mimc^- By 
repeating the same procedure to the remaining clusters, we obtain that max7r^ > max7r^, so that, if 
|II*| = m and p > p*, then < m. ■ 
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Following Theorem III.4 an optimal left-induced partition of cardinality (at most) m is also an optimal 
m-partition. Notice that for the computation of an optimal left-induced partition only the lengths p 
corresponding to the discontinuity points of IP need to be considered. Since each discontinuity point 
coincides with the distance between a pair of vertices, only n ^"~ 1 ^ values need to be tested. Therefore, an 
optimal left-induced partition can be computed with complexity 0(n 2 ). In what follows we describe an e- 
approximation algorithm with linear complexity for any e G M>o- Notice that e- approximation algorithms 
with linear complexity are often more convenient for a practical implementation than exact algorithms 
with higher complexity 0. 

We now present our algorithm for the computation of an optimal m-partition. Since the function p i-> 
|IP| is monotone and continuous, a bisection method is effective for finding its discontinuity points, and, 
therefore, for determining the shortest length of a left-induced partition of cardinality m. A bisection 
based procedure to compute an optimal left-induced partition is in Algorithm [2| where the function left- 
induced( {v i, . . . , v n },p) returns the left-induced partition defined in equation ([T]). We next characterize the 
convergence properties of Algorithm [2| 

Lemma IIL5 (Convergence of Algorithm [2]). Let G be a chain roadmap, and let Ii m denote an m- 
partition of G. Let p* = minn m dim(TL m ). Algorithm [2] with tolerance e returns a left-induced partition 
of dimension at most p* + e and cardinality at most m. Moreover, the time complexity of Algorithm^is 
O^log^" 1 )). 

Proof: Algorithm [2] searches for the minimum length p* that generates a left-induced partition of 



cardinality at most m. Because of Theorem III.4 , the length p* coincides with one of the discontinuity 
points of the function |IP|, and it holds p* E (0, 2v n /m). Indeed, p* > because m < n, and p* < 2v n /m, 
because (2v n /m)m > v n . Recall from ([2]) that |I1 P | > m for every p < p* , and that the function p \U P \ 
is monotone. Note that the interval [a, b], as updated in Algorithm [2j contains the value p* at every iteration. 
The length of the interval [a, b] is divided by 2 at each iteration, so that, after log 2 (f^), the value p* is 
computed with precision e. Since the computation of |n p | can be performed in 0(n) operations, the time 
complexity of Algorithm [2] is 0{n\og{e~ 1 )). 



As a consequence of Proposition III.2| and Theorem III.3, in what follows we only consider team 



trajectories whose image coincide with an m-partition. Therefore, for ease of notation, we use the set 
{Vx, . . . , V m } to denote both the image set of a team trajectory and an m-partition of the chain graph. 
We conclude this section with a summary of the presented results. 

Theorem IIL6 (Patrolling a chain graph at minimum refresh time). Let G be a chain graph with n 
viewpoints and let m < n be the number of robots. Let V be an optimal m-partition of G computed by 
means of Algorithm [2] with tolerance e. Let o? max be the dimension of V. A team trajectory with image 
V, and minimum refresh time 2o? max is computed as in Trajectory [7] Moreover, the time complexity of 
designing such trajectory is 0(n logfe -1 )). 



IV. Minimum refresh time and latency team trajectory on a chain roadmap 

The previous section considers the problem of designing team trajectories with optimal refresh time 
on a chain graph. In a patrolling mission it may be important for the robots to communicate with each 
other in order to gather information about the status of the entire environment. For instance, messages 
could be sent by a unit to ask for reinforcement, or to spread an alarm. Informally, we call latency of 
a team trajectory X, in short LT(X), the shortest time interval necessary for a message generated by 
any robot to reach all the other robots. In other words, given our communication model, the latency of 
a team trajectory is a measure of how fast a message spreads to all robots. In this section we describe 
team trajectories with minimum refresh time and latency. 

We now give a more formal definition of LT(X). Recall that, by assumption, two robots are allowed to 
communicate when they lie on two adjacent viewpoints. In a chain roadmap, for a message to reach every 
robot in the chain, every pair of adjacent robots needs to communicate. For i e {2, . . . , m}, let $j denote 
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the union of the set of times at which the robots i — 1 and i communicate and {0}. The up-latency of X, 
in short LT up (X), is the longest time interval between any two consecutive communications between the 
robots 1, 2 and m — 1, m. Precisely, 

LT up (X) = max min t m — i 2 , 

t2G*2 t m e# m (t 2 ) 

where $ m (t 2 ) = {t m G $ m | 3t 3 G $ 3 , . . . ,t m _i G $ TO _i,t 2 < t 3 < • • • < t m } U {!>}. Analogously, we 
call down-latency the quantity 

LT down (X) = max min t 2 -t m , 

where $ 2 (t m ) = G $2 | 3t 3 G $3, . . . ,t m -i e $m-i,t 2 > *3 > • ■ • > W U {T f }. Finally, we define 
the latency of a team trajectory as 

LT(X) = max{LT up (X),LT down (X)}. 

Notice that our definitions of latency hold for m > 2, and that, if m — 2, then we have LT up (X) = 
LT down (X) = LT(X) = for every team trajectory X. We envision that the up- and down-latency 
performance criteria should be adopted when it is of interest to report the status of the monitored area 
to a base station located at one end of the chain environment. The latency minimization problem is more 
appropriate for fully distributed scenarios. In this section we design synchronized team trajectories with the 
following two features. First, since a minimum refresh time trajectory is determined by an optimalpartition 
of the chain graph, we aim at finding team trajectories associated with the same optimal partition^ Second, 
we design synchronized team trajectories with minimum up-latency (resp. down-latency) or latency. 



A. Lower bound and optimal team trajectory for up-latency 

We start by showing a lower bound for LT up (X) and LT down (X). Recall that, for a partition {Vi, . . . , V m }, 
we have k = mxn veVi v, r* = max,,^ v, d { = r { - l u and d max = max ie{lj ... im} d { . 

Lemma IV.l (Up-latency lower bound). Let G be a chain roadmap, and let {Vi, . . . ,V m } be an m- 
partition of G. The latency of a team trajectory with image {Vi, . . . , V m } is lower bounded by Y1T=2 di- 

Proof: The proposition follows from the fact that the robots speed is bounded by 1, and that the 
robots need to travel their segment to communicate with the neighboring robots. ■ 



For the up-latency of a team trajectory to equal the lower bound in Lemma IV.l each robot i needs 
to transfer a message from robot % — 1 to robot % + 1 in time dj. In order to do so, each robot i needs to 
communicate with its neighbor i + 1 as soon as Xi(t) = 

Theorem IV.2 (Patrolling a chain graph at minimum refresh time and minimum up-latency). Let G be a 

chain graph with n viewpoints and let m <n be the number of robots. Let V be an optimal m-partition 
of G computed by means of Algorithm [2] with tolerance e. Let d max be the dimension of V, and let di be 
the length of the i-th cluster. A team trajectory with image V, minimum refresh time 2<i max , and minimum 
up-latency ^^=2 4/ zs computed as in Trajectory 3 Moreover, the time complexity of designing such 
trajectory is 0(n log(e -1 )). 

Proof: The theorem follows by observing that the trajectory is 2d msx -periodic, and that no robot i 
waits at to communicate with the neighboring robot i + 1. The up-latency equals the lower bound 



in Lemma IV. 1 , and it is therefore minimum. Regarding the computational complexity, notice that it 



is determined by the computation of the optimal m-partition V, and hence, by Lemma III.5 it equals 
0(nbg(e" 1 )). ' ■ 
An example of a team trajectory with minimum refresh time and minimum up-latency is in Fig. |4j 
Finally, observe that the minimization of the down-latency can be achieved in an analogous way. 
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Trajectory 3: Minimum base-latency team trajectory (i-th robot) 

Input : U := min^gVi v, Vi := max„ 6 y, v, di := Ti — U, d max := max., rj — lj\ 
Require : optimal partition of the chain graph; 

Set : to := - Y^i d h k £ N h-> T<(fc) := 2fcd raiix + t + E*-=i dj! 

i: := i< for t <t< 2i(0) and for Ti(fc) + 2d; < £ < Ti(fc + 1); 

2: Ki(£) := n for t := Ti(ft) + d*; 
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Fig. 4. A team trajectory with minimum up-latency generated with the procedure in Trajectory [3] Notice that each robot i communicates 
with the neighboring robot i + 1 as soon as Xi(t) = A team trajectory with minimum down-latency is obtained by shifting in time the 
trajectory of robot 1, in a way that robot 2 communicates with robot 1 as soon as x%(t) = h- 



Remark 3 (Simplifications for clusters of equal length). Let {Vi, . . . ,V m } be an optimal m-partition, 
and suppose that <i max = r« — U for all i £ {1, . . . , m}. Then the up-latency and the down-latency can be 
made minimum and equal to (m — 2)c/ max by arranging the robots trajectories to be in opposite phase. 
Specifically, for k = 0, 2, ... , we set 

Xi(2kd max ) = Zj, Xj((2fc + l)<i m ax) = f^, 

if i is odd, and 

Xi(2kd max ) Tj, Xj((2A; -)- l)<i max ) Zj, 



if i is even. Because of Lemma IV. 1 the above trajectory has minimum latency. This particular case was 
studied in /|2]/. 



B. Lower bound for latency 

We now consider the minimization of the latency criterion, and we restrict our attention to periodic 
team trajectories. To be more precise, let {Vi, . . . , V m } be an optimal m-partition of the environment, and 
let d max denote the longest length of the clusters. We aim at finding a 2<i max -periodic team trajectory with 
image {Vi, . . . ,V m } and minimum latency. Notice that, by imposing a periodicity of 2d max , the refresh 
time of the trajectory, if finite, is also minimized. 

We start by considering the pedagogical situation in which di + d i+ % > d max for allz£{l,...,m — 1}. 
In the next Lemma, we show that the frequency of message exchange among the robots is limited by the 
periodicity of the trajectory. 

Lemma IV.3 (Frequency of message exchange). Consider a 2d max -periodic team trajectory, and let di + 
d i+1 > d max for all i £ {1, . . . , m — 1}. For any t £ [0, T f — 2d max ] and for any i £ {2, . . . , m — 2}, there 

'We focus on this family of trajectories, and we leave the more general optimization problem as the subject of future research. However, 
this family of trajectories is believed to contain an (unconstrained) optimal solution. 
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Fig. 5. As stated in Lemma [lV.3| at most one communication sequence is possible within each period of length 2d raax . Here di > d mla . 

exist no two distinct sequences ti, t i+ \, t i+2 and U, with tj, tj e j — i, i + 1, % + 2, such that 

t < ti < < t i+2 < t + 2d max 
t <ti < t i+ i < t i+2 < £ + 2<i max 

t+1 < t 

Moreover, for any i G {m, . . . , 4}, ?/zere exxsf no two sequences ti, £j_ 2 and tj, £j_ 2 ?j £ 

$ 3 , j = i,i — l,i — 2, such that 



max 



t<U< U-i < ti_2 < t + 2d ; 
t < ti < ti-i < < £ + 2<i max 

Proo/- Since + d i+2 > d m ax, it follows max{d i+ i,(l 1+2 } > d max /2. Let > d max /2. By 
contradiction, if two distinct sequences U,t i+ i,ti + 2 and U,ti + i,t i+ 2 exist, with t < U < t i+ i < U + 2 < 
t + 2d max , t <ii < < i i+2 < t + 2d max , t i+1 < i{ and t i+2 < then the (i + l)-th robot travels its 
cluster four times. Since the speed of the robots is bounded by one, robot i + l cannot travel its cluster 
four times in a period of 2d max (cf. Fig. [5]). The second part of the theorem follows from an analogous 
reasoning. ■ 
Notice that in the above Lemma the index i belongs to the set {2, . . . ,m — 2} (resp. {m, . . . ,4}) 
because we consider 3 consecutive communication instants, and because $j denotes the sequence of times 



at which robots % — 1 and % communicate. Because of Lemma IV.3 in a 2<i max -periodic team trajectory 
with di + di + i > d max only one message can be passed from robot i to robot i + 3 in a period of time 
of 2d max . This limitation determines a lower bound on the latency of a periodic trajectory. Notice that 
eventual communication instants tj E $j, with t < tj < ti < t i+ i, do not influence the latency, since all 
information can be transmitted at time U without affecting the latency performance. 

Lemma IV.4 (Latency lower bound, simple case). Let X be a 2d max -periodic team trajectory with di + 
d i+1 > d max for all i G {1, . . . , m — 1}. Then 

minLT(X) > (m — 2)d max . 



Proof: Because of Lemma IV.3, a message can be transfered from robot i to robot i + 3 in at most 
2e?max instants of time, by traveling the clusters V i+ \ and V i+2 . Without losing generality, we let d max be 
the time to pass a message form % + 1 to % + 2, and from i + 2 to i + 3. Notice that the same reasoning 
holds also for the time to transfer a message from i + 3 to i. Therefore, the latency is lower bounded by 
(m - 2)d max . M 
We now consider the situation in which an optimal m -partition does not verify the constraint di + d i+ i > 
dmax for allie{l,...,m — 1}. Intuitively, for what concerns the latency, two consecutive clusters with 
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length di and may be regarded as one single cluster of length d{ + d i+ \. Therefore, in order to use 
Lemma IV.3 and Lemma IV.4 we partition the clusters {Vi, . . . , V m } into groups such that the sum of the 
length of the clusters of two consecutive groups is greater than d max . In other words, let V r = {f 1 , . . . , f m } 
be the set of right-extreme viewpoints of the partition {Vi, . . . ,V m } defined recursively as 

r := h, 

j 

max {rj V4< d max }, i = 1, 

je{l,...,m} 



,m, 



k=pi 



where pi = min{l, . . . , m} such that l Vi > jvi. Let Vi — {h, . . . , l m } be the set of left-extreme viewpoints 
defined recursively as 



h '■— h, 

U := min {lj \ lj > fj_i}, i — 2, . . . , m. 
je{i,...,m} 

Additionally, define the set of aggregated clusters V = {V±, . . . , V m }, where V contains all the clusters 
within li and fj, and let di be the sum of the length of the clusters in Vi. 

Lemma IV.5 (Latency lower bound, general case). Let X be a 2d max -periodic team trajectory with image 
{Vi, . . . , V m }, and let m be the number of aggregated clusters. Then, 

minLT(X) > (m - 2)d max + {d\ - d x ) + (d m - d m ). 

Proof: Consider the clusters defined by the right extreme viewpoints, and notice that they verify 
di + di + \ > d max . Then, the Theorem follows from Lemma IV.4, and from the fact that the minimum 
latency on the image {Vi, . . . , V m } equals the minimum latency on the image {Vi, . . . , V m }. The terms 
d\ — d\ and d m — d m are due to the fact that we are interested in delivering a message from robot 1 to 
robot m in the original configuration, and not on the aggregated chain. ■ 



C. Optimal team trajectory for latency 

A team Trajectory with minimum refresh time and minimum latency is formally presented in Trajectory 
|4} where we specify the instants of time at which a robot changes its velocity, and we assume that it 
moves at maximum speed otherwise. An example is reported in Fig. [6[ and here we give an informal 
description. Let {Vi, . . . , V m } be an optimal m-partition of a chain graph, and let {V\, . . . , V m } be the set 
of aggregated clusters. Recall that Vi is a subset of {Vi, . . . , V m }, and that the sum of the length of two 
consecutive elements is larger than d m3X . The procedure in Trajectory [4] (lines 6 — 11 and 18 — 23) is such 
that the robots in the same group behave as a single robot assigned to the whole set of viewpoints. In other 
words, the motion of the robots in the same group is determined by a token passing mechanism, in which 
robot 2 + 1 moves towards r i+ % only when Xi{t) = and Xi(t~) ^ r i5 and, analogously, robot i moves 
towards k only when Xi(t) = k and Xi(t~) ^ U. Instead, lines 1 — 5 and 12 — 17 in Trajectory [4] guarantee 
the transfer of one message in a period of 2<i max between three consecutive groups, and, consequently, the 
minimization of the latency. Indeed, since the sum of the length of two consecutive groups is larger than 
dmax, because of Lemma |IV.3[ no more than one message can be transferred between three consecutive 
groups in a period of 2<i max . 

Theorem IV.6 (Patrolling a chain graph at minimum refresh time and minimum latency). Let G be a 

chain graph with n viewpoints and let m < n be the number of robots. Let V be an optimal m-partition 
of G computed by means of Algorithm [2] with tolerance e. Let d max be the dimension of V, and let d\ 
(resp. d m ) be the length of the first (resp. last) cluster in V. Let m be the number of aggregated clusters, 
and let d\ (resp. d m ) be the length of the first (resp. last) aggregated cluster. A team trajectory with image 
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Trajectory 4: Minimum refresh time and latency trajectory (i-th robot) 

Input : U := min^g^ v, r$ := max^g^ v, d max := maxj Vj — lj, R := set containing the identifiers of the robots in the same 

group as i, k := index of the aggregated cluster {1, . . . , m}; 
Require: optimal partition of the chain graph; 

case U. is a left-extreme 
if k is odd then 
| Xi(t) = h with t = 0, 2d max , 4d max , . . . ; 

else if k is even then 
| Xi(t) = h with t = d max , 3d m ax, 5d max , . . . ; 



case U is not a left-extreme 

fit : = £ieit,i<« 
if fe w odd then 

|_ = for r — 5j < t < r + Si, with r = 0, 2d max , 4d max , . . . ; 

else if k is even then 
|_ Xi(t) = /j for r - <5; < t < r + 5,, with r = d max , 3d milx , 5d raax , . . 

case ri is a right-extreme 

Si '.— dmax XyjgJ? 

if fe is odd then 

| a;»(t) = r l for t - Si <t < r + 5i, with r = dmax, 3d max , 5d max , . . 

else if k is even then 
| Xi(t) = n for r — Si < t < t + Su with r = 0, 2d max , 4d max , . . . ; 

case ri is not a right-extreme 
Si :— d max — J2jeRj<i ^J' 
if is odd then 

|_ Xi(t) = n for t - Si < t < t + Si, with r = dmax, 3d max , 5d max , . . 

else if k is even then 
|_ Xi(t) = n for r — Si < t < r + Si, with r = 0, 2d max , 4d max , . . . ; 



V, minimum refresh time 2d max , and minimum latency (m — 2)<i max + d\ — d\ + dm — d m is computed as 
in Trajectory^ Moreover, the time complexity of designing such trajectory is 0(n log(£: _1 )). 

Proof: By inspection, the team Trajectory described in Trajectory [4] is 2<i max -periodic, and therefore 
it has minimum refresh time. Moreover, by construction, the communications at the extreme viewpoints 
happen every 2d max instants of time, so that the latency is equal to (rh — 2)d max + d\ — d\ + dm — d m , and 
hence, by Lemma IV.5, minimum. Regarding the computational compl exity, notice that it is determined 
by the computation of the optimal m-partition V, and hence, by Lemma HI. 5 it equals 0(n log(e -1 )). ■ 



V. Distributed synchronization algorithm on a chain roadmap 

In the previous sections we have shown that, for the computation of a minimum refresh time and latency 
team trajectory, first an optimal m-partition of the roadmap needs to be found, and, then, a synchronization 
of the motion of the robots needs to be achieved to ensure communication between neighboring robots. 
The distributed computation of an optimal m-partition follows directly from Algorithm [2[ by letting the 
robots compute the left-induced partition of length I in a distributed way. A simple solution consists of 
the following three operation: 

(i) the robots gather at the leftmost viewpoint, and 

(ii) determine the cardinality of the team and elect a leader; 

(iii) the leader computes an optimal left-induced partition, and assigns a different cluster to each robot. 
Notice that, by only assuming the capability of detecting the viewpoints in the roadmap (in particular 
the leftmost and rightmost viewpoint) the robots can distributively compute an optimal partition. Indeed, 
the leader can simply travel the roadmap, and evaluate if, for a given length p, the cardinality of the 
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viewpoints 




Fig. 6. A team trajectory with minimum refresh time (2d max ) and minimum latency ((m — 2)<f max + d\ + di — d\ + d% — d$ = 
di) is here reported. A message with delivery time equal to di is reported in red. The right-extreme viewpoints are {fi,f2}, the left- 
extreme viewpoints are {h,li}, and rh — 2. The represented team trajectory has minimum refresh time (2d raax ) and minimum latency 
((m — 2)cf max + (di + di) — di + ds — d3 = di). A message with delivery time equal to di is reported in red. 



corresponding left-induced partition is greater, equal, or smaller than the cardinality of the team. We 
believe that, by simple programming operations, the above procedure can be improved, so as to handle 
reconfigurations due to addition or removal of robots in the team. In this section we focus instead on 
the design of a distributed feedback law to synchronize the motion of the robots so as to guarantee the 
minimization of the latency of the trajectory. 

Recall that Xi(t) denotes the position on the chain of the robot i at time t. Moreover, let the i-th 
cluster of an optimal m-partition be delimited by k and ?v Let dir.; 6 { — 1,0, 1} denote the direction 
of motion of robot i. Let c-time denote the current time of simulation, let a-time be the time at which 
a robot arrives at his right extreme. Let n-meet(ij) be a function that returns the number of times that 
the robots i and j have communicated. Let Timer(S) be a function that returns 1 after a time interval of 
length 5. An algorithm for the robots to distributively converge to a minimum refresh time and latency 
team trajectory is in Algorithm [5] It should be noticed that the algorithm assumes the knowledge of an 
optimal partitioning of the chain graph, and of the left- and right-extreme sets. 

Algorithm [5] is informally described as follows. First, the velocity of a robot changes only when its 
position coincides with the leftmost or the rightmost of the assigned viewpoints. When a robot reaches 
an extreme viewpoint, it waits until a communication with the neighboring robot happens (lines 7 — 8). 
This mechanism determines the feedback behavior of our procedure. The behavior of a robot after a 
communication is determined by the lines 9 — 21, which reproduce the optimal behavior described in 
Trajectory |4} To be more precise, the function Token(i,j) coordinates the motion of the robots in the same 



group (see Section IV-C), so that they move as if there was a single robot sweeping the viewpoints in 
the same group. The function Timer(Si), instead, ensures the maximum frequency of messages exchange 
between consecutive groups, so as to minimize the latency of the resulting trajectory. 

Theorem V.l (Optimal team trajectory). Let X be the team trajectory generated by Algorithm [5] There 
exists a finite time after which X has minimum refresh time and latency. 

Proof: Let {Vi, . . . , V m } be an optimal m-partition, let V r — {fi , . . . , f c * } be the set of right-extreme 
viewpoints, and let Vi — {It, . . . , lc*} be the set of left-extreme viewpoints. Let R denote the set of robots 
patrolling a viewpoint between ^ and fj, where k and fj are as previously defined. First, notice that robot 
1 (resp. m) sets xi(t) = 1 (resp. x m (t) = —1) as soon as xi(t) = h (resp. x m (t) = r m ). Second, the 
function Token(i,j) guarantees that, when i,jeR communicate, exactly one robot among i,j maintains 
a zero velocity and in an alternate way. Therefore, after a finite time T i? independent of the initial robots 
motion direction, the velocities of the robots in R are such that, upon communication, Xi(t) = 
and ij-i(t) = 0. In other words, after Tj, the robots in R behave as a single robot sweeping the segments 
between ^ and fj. Finally, the function Timer(r) and the parameter r guarantee that the communications 
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Algorithm 5: Minimum refresh time and latency team trajectory (i-th robot) 



v, cL 



Input : U := min„ g Vi v, ri := max^gv; 
group as i; 

Set : Si := (d raax - £\ eB d lV 2 ' dir ' G i 1 ' 

Require : optimal partition of the chain graph; U < Xi(0) < r 

while true do 

case i = 1 and x±(t) = l\ 
[ diri := 1; 

case i = m and x m (t) = r m 
|_ dir m := -1; 

case (Xi(t) = ?j and Xi—i(t) 7^ rj_i) or 
|_ diu := 0; 

case = and = rj_i 

if G V; then 

receive r from robot i — 1; 
dki := Timer(r); 

else if /i ^ V; then 

|^ dki := Token(i, i ~ 1); 

case Xi(t) — Ti and Xi+\{t~) — U+i 

r := max{0, a-time + Si — c-time}; 
send t to robot i + 1; 
if ri G t4 then 
I di^ := — Timer(5i + r); 

else if ri V T then 

J dki := Token(i, i + 1); 

±i{t) := dm; 



maxj rj — lj, R := set containing the identifiers of the robots in the same 



(Xi(t) = n and x i+ i(t) 7^ h+i) 



Function Token(ij) 

i: case i > j if n-meet(ij) is even then return 1 
2: case i < j if n-meet(ij) is odd then return -1 
3: otherwise return 



at the extreme viewpoints happen every 2d max instants of time. We conclude that the trajectory generated 
by Algorithm [5] converges in finite time to a team trajectory with minimum refresh time and latency. ■ 
It should be observed that, differently from the team trajectories presented in the previous sections, 
Algorithm [5] contains a feedback procedure to synchronize the motion of the robots. Our algorithm is 
independent of the initial configuration, and, as it is shown in the next section, it is robust to a certain 
class of robot failures. 

VI. A CASE STUDY 

Given the theoretical nature of this work, in this section we propose a simulation study to show the 
effectiveness of our procedures. For our simulations, we use the Matlab® simulation environment, and 
we model the robots as holonomic vehicles of zero dimension. The communication edges and the motion 
paths are described by a given roadmap. 

Consider the chain roadmap with 30 viewpoints in Fig. |7J Suppose that 10 robots are assigned to the 
patrolling task. In order to obtain a team trajectory with minimum refresh time, an optimal 10-partition 
of the roadmap is computed (cf. Fig. [7]). Additionally, to obtain a minimum latency team trajectory, each 
robot is assigned to a different cluster, its initial position is chosen randomly inside its cluster, and its 
velocity is initialized randomly. The motion of each robot is then determined by Algorithm [5] The resulting 
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Fig. 7. The figure shows a robotic roadmap with sensor coverage and communication connectivity of a two floors building. Crossing edges 
corresponds to corridors at different floors. A 10-partition of the roadmap is here reported. The dashed edges are not traveled by any robot. 
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Time (s) 

Fig. 8. For the roadmap of Fig. [7] the team trajectory obtained with Algorithm [5] is reported here. The dashed lines denote a communication 
among two neighboring robots. At time t = 200 the robots have synchronized their motion, and from that moment up to time t — 300 
the team trajectory has minimum refresh time and latency. From time t = 300 up to time t — 400, robot 7 undergoes a temporary failure, 
causing all the other robots to lose synchronization. The team of robots synchronizes again when robot 7 resume its motion. At time t = 620, 
the failure of robot 7 is detected by the remaining robots, which compute and synchronize on a new partition. 

team trajectory is in Fig. [8j where the team of robots synchronize on a minimum refresh time and latency 
team trajectory after a finite transient. 

We now test the robustness of our synchronization procedure. As first case study, we consider a 
temporary stopping failure, in which a robot stops for a certain time interval. For instance, suppose 
that robot 7 stops from time 300 up to time 400 (cf. Fig. [8]). Notice that, after the failure, each robot j, 
with j < 7, gathers at rj, and each robot k, with k > 7, gathers at l k waiting for a communication with 
the corresponding neighboring robot. As soon as robot 7 resumes its functionalities, the team of robots 
recover the desired synchronization. Notice that the transient failure of robot 7 can be easily detected by 
its neighbors by means of a timer mechanism with a predefined threshold. 

As a second case study, we let the robots actuation be affected by noise, so that the speed of the robots 
becomes a random variable with a certain distribution^] Precisely, let ±i = dir^ + Wi be the equation 
describing the dynamics of robot i, where Wi is a zero mean Gaussian variable with variance a 2 [(m/s) 2 ]. 
We let o 2 E {0, 0.02, . . . , 0.5} and we run 100 simulations for each possible value of a 2 on the roadmap 
of Fig. [7] The refresh time and the latency of the team trajectories obtained with Algorithm [5] are plotted in 



2 The case in which a robot fails at seeing a neighboring robot for a certain interval of time can be modeled analogously. 
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Fig. 9. For the roadmap in Fig. [7] 100 team trajectories have been generated with Algorithm [5] for each value of a 2 . The average refresh 
time and the average latency are plotted as continuous lines in Fig. |9(a)| and Fig. |9(b)| respectively. The bars indicate the minimum and the 
maximum value of the refresh time and the latency, respectively. Notice that the average team trajectory performance degrade gracefully 
with the noise variance. 



Fig. 9(a) and in Fig. 9(b) respectively, as a function of a 2 . Note that the performance degrade gracefully 
with the noise magnitude. 

As third and final case study, we consider the situation in which a robot definitively stopsj^] The 
remaining robots need to compute a new optimal partition and to synchronize in order to guarantee an 
optimal patrolling of the environment. Notice that for the computation of such a partition by Algorithm 
[2] the chain graph and the number of the robots is required. Suppose that the failure of the robot 7 is 
detected at time 620 by the well-behaving robots, and assume that each robot knows the chain roadmap 
and the number of operative robots. Algorithm [2] and Algorithm [5] allow the team to synchronize on a 
new team trajectory with minimum refresh time and latency. Notice that the initial and the final partitions 
do not coincide. 



VII. Approximation algorithms and heuristics for general roadmaps 

The problem of designing minimum refresh time and latency team trajectories on a chain roadmap 
has been discussed. In this section we consider the more general cases of tree and cyclic roadmap, 
we characterize the computational complexity of determining optimal trajectories, and we describe two 
approximation methods with performance guarantees. The results we are going to present are intended 
for a team of more than one robot. Indeed, if only one robot is assigned to the patrolling task, then a 
minimum refresh time trajectory follows from the computation of the shortest tour through the viewpoints, 
for which efficient approximation algorithms already exist ll2~8l . 



A. Minimum refresh time team trajectory on a tree roadmap 

Let T = (V, E) denote an undirected, connected, and acyclic roadmap (tree). Recall that a vertex path 
is a sequence of vertices such that any pair of consecutive vertices in the sequence are adjacent. A tour 
is a vertex path in which the start and end vertices coincide, and in which every vertex of T appears at 
least once in the sequence. A depth-first tour of T is a tour that visits the vertices V in a depth-first order 
Il2~9~ll . Let DFT(T) denote the length of a depth first tour of T. Notice that the length of a depth-first tour 
of a connected tree equals twice the sum of the length of the edges of the tree, and that any depth-first 
tour is a shortest tour visiting all the vertices. We now show that, for the case of tree roadmap, the set 
of cyclic-based and partition-based trajectories described in [fT51 does not contain, in general, a minimum 
refresh time trajectory. Recall that in a cyclic -based strategy the robots travel at maximum speed and 
equally spaced along a minimum length tour visiting all the viewpoints. Consider the tree roadmap of 
Fig. |10(a)[ and suppose that two robots are assigned to the patrolling task. Clearly, the minimum refresh 
time is 2s, while the refresh time of a cyclic strategy equals l+e. Consider now the tree roadmap in Fig. 



10(b) where the edges have unit length, and assume that two robots are in charge of the patrolling task. 



3 The case of additional robots joining the team is handled in a similar way. 
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Fig. 10. Two examples of tree roadmap. For the roadmap i n Fig. 1 10(a)| the family of cyclic based trajectories does not contain a minimum 
refresh time team trajectory. Instead, for the roadmap in Fig. |10(b)| the~family of partitioned based trajectories does not contain a minimum 
refresh time team trajectory. 



Observe that any partition of cardinality 2 contains a chain of length 2, so that, since only one robot is 
assigned to each cluster, the minimum refresh time that can be obtained is 4. Suppose, instead, that the 
robots visit the vertices of the roadmap as specified in Table [j] where x(t) denotes the position of a robot 
at time t. Since the refresh time of the proposed trajectory is 3, we conclude that neither the cyclic-based 
nor the partition-based strategy may lead to a minimum refresh time team trajectory on a tree roadmap. 



TABLE I 



Robot 


x(0) 


x(l) 


x(2) 


x(3) 


x(4) 


x(5) 


x(6) as (7) 


1 


Vl 


V2 


«4 


V2 


vz 


«2 


Vl 


2 


V2 


«3 


Vl 


Vl 


V2 


Va 


V2 



We now introduce some definitions. Let X be a team trajectory on the tree roadmap T. We say that 
the edge (vj,v z ) E E is used by X if there exists i E {1, . . . , m} and (ti,t 2 ) E [0,RT(X)] 2 such that 
Xi(ti) = Vj and Xi(t 2 ) = v z , and it is unused otherwise. Note that, because in a tree there exists only 
one path connecting two vertices, the above condition ensures that the edge (J, z) is traveled by the robot 
i. Let E denote the set of unused edges, and let Ft be the forest obtained from T by removing the 
edges E from E, i.e., the collection of vertex-disjoint subtrees {Ti, . . . , T^}, with Tj = (Vi,Ei), such 
that V = uf =1 Vi and Ei C E, for each i E {1, . . . , k}. Let rrii be the number of robots that visit at 
least one vertex of Tj in the interval [0, RT(X)], and note that > in a finite refresh time trajectory. 
Let M = {mi, . . . ,m k }. Notice that the same subtree collection can be associated with different team 
trajectories. We say that a team trajectory is efficient if its refresh time is the smallest among all the team 
trajectories associated with the same subtree collection. 

Theorem VII. 1 (Efficient team trajectory). Let (F T ,M) be the subtree collection associated with the 
team trajectory X on the tree roadmap T, where F T = {T 1; . . . ,T k }, and M = {mi, . . . , m^}. Then, X 
is efficient if 

RT{X) = max DFT (T^/rrij. 
je{i,...,k} 

Proof: Let i E {1, . . . , k}, and let m ; be the number of robots assigned to Tj. Notice that the robots 
in Tj travel, in total, at least DFT(Tj) to visit all the vertices. Since the speed of the robots is bounded 
by 1, the smallest refresh time for the vertices of Tj is DFT (T^/m*. ■ 
An efficient team trajectory, can be computed with the following procedure. See Table [I] for an example. 

Lemma VII.2 (Efficient team trajectory computation). Let (Ft, M) be a subtree collection of a tree 
roadmap, where Ft = {Xi, . . . , T^}, and M = {mi, . . . , m k }. An efficient team trajectory is as follows: 
for each i E {1, . . . , k}, 

(i) compute a depth-first tour n of T it 

(ii) equally space robots along n, and 
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(iii) move the robots clockwise at maximum speed on r». 

Proof: Since every vertex of Tj G i<V appears at least once in a depth first tour r< of Tj, and the 
robots move with maximum speed and equally spaced along n, every vertex is visited at most every 
DFT(Tj)/mj. The statement follows. ■ 
Let P[m) be the partition set of m, i.e., the set of all the sequences of integers whose sum is m. The 
following problem is useful to characterize the complexity of designing minimum refresh time trajectories 
on a tree roadmap. 

Problem 2 (Optimal subtree collection). Let T be a tree roadmap and let m be the number of robots. 
Find a subtree collection (Ft, M) that minimizes maxj6{i ) ... ) |ir T |} DFT(Tj)/rrij subject to M e P(m) and 
\F T \ = \M\. 

Lemma VII.3 (Equivalent problem). For the case of a tree roadmap, the Team refresh time problem and 
the Optimal subtree collection problem are equivalent. 



Proof: As a consequence of Theorem VII. 1 , the minimum refresh time on a tree roadmap T can be 
written as mm^ FT:M ) maXj e /i jm. DFT (Tj)/rrij, where (F,M) is a subtree collection of T, and \M\ = 
\Ft\ — k < m. It follows that a solution to Problem [T] can be derived in polynomial time from a solution 



to Problem [2] by using the procedure described in Lemma VII. 2 Suppose now we have a solution to 



Problem [TJ Then an optimal subtree collection follows from the identification of the unused edges. We 
conclude that the two optimization problems are equivalent. ■ 
We now state our main result on the design of minimum refresh time team trajectory on a tree roadmap. 

Theorem VII.4 (Computing a minimum refresh time team trajectory on a tree). Let T be a tree roadmap 
with n vertices, and let m be the number of robots. A minimum refresh time team trajectory on T can be 
computed in 0((m — l)!n) time. 

Proof: Recall from ll30l that an optimal subtree collection can be computed in 0((m — l)!n). Then, 
by using Lemma | VII.3 1 and Lemma VII2[ the claimed statement follows. ■ 



As a consequence of Theorem VI 



.4 the problem of designing minimum refresh time team trajectories 
on a tree roadmap is computationally easy for any finite number of robots. In our design procedure, we 
first compute an optimal subtree collection of the given tree, and then we schedule the robots trajectory 
according to Lemma fVTI.2 



B. Minimum refresh time team trajectory on a cyclic roadmap 

In this section we propose two approximation methods for the Team refresh time problem in the case 
of a cyclic, i.e., not acyclic, roadmap. These solutions are obtained from a transformation of the cyclic 
roadmap into an acyclic roadmap. 

Let G = (V, E), with \V\ = n, be an undirected and connected roadmap. Note that there exists an open 
tour r with at most 2n — 4 edges that visits all the vertices^] We construct a chain roadmap F from r by 
doubling its repeated vertices and edges, so that T has at most 2n — 3 vertices and at most 2n — 4 edges, 



and such that the length of the i-th edge of F equals the length of the i-th edge of r (cf. Fig. 11 ). Our 
first approximation method consists of applying Algorithm [5] to an optimal m-partition of F. 

Theorem VII.5 (Performance ratio). Let G be a connected roadmap, let n be the number of vertices of 
G, and let 7 be ratio of the longest to the shortest length of the edges of G. Let i?7* be the minimum 
refresh time on G. Let t be an open tour with 2n — 4 edges that visits all the n vertices, and let F be 
the chain roadmap associated with r. Let RT^ be the minimum refresh time on F. Then 

R7*r < \ ~~) 8 7^ r - 



4 An open tour with at most 2n — 4 edges that visits all the vertices can be constructed starting from a leaf of a spanning tree of G. 
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Fig. 11. The figure shows a chain roadmap associated with the cyclic roadmap of Fig. [T] Notice that the cycles in Fig. [T] have been broken. 
Moreover, 3 vertices (3 edges) of Fig. [T] are repeated twice in the chain. 




£ 1 1 1 1 1 
• • • • • • 



1 e e 1 1 1 



Fig. 12. A tree roadmap and two corresponding chain roadmaps. If the number of robots is 4, then the performance ratio RTf/RT* grows 
with e _1 . 



Proof: Let w be the shortest length of the edges of G, and note that the length of T is upper bounded 
by 2{n — 2)^w. It follows that RT r < 4< "^f 7 ~ - Since m < n by assumption, some robots need to move 
along G for all the viewpoints to be visited. Because each robot can visit only a vertex at a time, at least 
| — — l] steps are needed to visit all the vertices of G, and therefore RT* > | — — l] w > \-w_- By 
taking the ratio of the two quantities the statement follows. ■ 
It should be noticed that, when 7 grows, the performance of our procedure might degrade. For instance, 



suppose that the roadmap is as in Fig. 12 and suppose that 4 robots are assigned to the patrolling task. 
As long as e < 1, a minimum refresh time strategy requires one robot to patrol the viewpoints {t>i,t>2}, 
while the second, third, and fourth robot stay on the viewpoints t> 3 , t> 4 , and t> 5 respectively. It follows that 
RT* = 2e. On the other hand, an optimal m-partition of any chain graph associated with a tour that visits 
all the viewpoints has dimension at least 1. Consequently, the refresh time of the team trajectory obtained 
with Algorithm [5] equals 2, and the ratio RTp/RT* grows proportionally to e _1 . 

We next describe a polynomial time constant factor approximation algorithm for the minimum refresh 
time problem. Given a roadmap G = (V, E) and a positive integer k < \V\, we define a path cover of 
cardinality k as the collection of paths {pi, . . . ,pk} such that V C {J* =1 Pi. Let the cost of a path equal 
the sum of the lengths of its edges. The min-max path cover problem asks for a minimum cost path 
cover for the input graph, where the cost of a cover equals the maximum cost of a path in the cover. The 
following result is known. 

Theorem VII.6 (Min-max path cover 0). There exists a 4- approximation polynomial algorithm for the 
NP-hard min-max path cover problem. 



Following Theorem |VIL6 , given a graph G, there exists a polynomial time algorithm that computes a 
path cover of G with cost at most 4 times greater than the cost of any path cover of G. We now state 
our approximation result for the NP-hard Team refresh time problem. 



Lemma VIL7 (8-approximation refresh time). There exists an 8- approximation polynomial algorithm for 
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Fig. 13. The picture shows the trajectories of 7 robots for the patrolling of a part of the UCSB campus. The viewpoints (red circles) are 
chosen so as to provide sensor coverage of the whole area. For the design of the patrolling trajectory, a roadmap is first constructed as 
described in Section [n] Then, a path covering of cardinality 7 is computed with the procedure in (7), and each robot is assigned a different 
path. Finally, the trajectory of each robot consists in sweeping at maximum speed the tour obtained by shortcutting the assigned path. The 
refresh time of the proposed team trajectory is proven to be within a factor of 8 of the minimum refresh time for the given roadmap. 



the NP-hard Team refresh time problem. 

Proof: Let {pi, . . . ,p m } be a 4-approximation path cover of the graph G. Note that the length of 
each path is within 4RT*. Indeed, in a minimum refresh time team trajectory starting at time and with 
unitary velocity, every vertex is visited within time RT*. Let X be a team trajectory obtained by letting 



robot i sweep at maximum speed the path p^. Clearly, RT(X) < 8RT*. Because of Theorem VII.6 the 
team trajectory X can be computed in polynomial time. ■ 



Following Lemma VH7[ for any given roadmap and any number of robots, a team trajectory with 



refresh time within a factor of 8 of the optimal refresh time can be constructed by computing a path 
covering of the roadmap, and by assigning a different path to each robot. An example is in Fig. [T3j a 
movie of which is included in the multimedia material. 

Remark 4 (Improving the team trajectory). Several existing heuristics can be used to improve upon the 
trajectories in Fig. [O] For instance, since the robots move in a metric space, shortcutting techniques may 
be applied 071/ . Because these heuristics do not guarantee an improvement of the optimality gap of our 
trajectories, they are not considered in this work, and they are left as the subject of future investigation. 



VIII. Conclusion and future work 

The design of team trajectories to cooperatively patrol an environment has been discussed. After defining 
the problem and the performance criteria, we analyze the computational complexity of the design problem 
as a function of the shape of the environment to be patrolled. For the case of a chain environment, we 
describe a polynomial algorithm to compute minimum refresh time and latency team trajectories. For the 
case of a tree environment, under the technical assumption of a constant number of robots, we identify 
a polynomial time algorithm to compute minimum refresh time team trajectories. Finally, the general 
case of cyclic environment is shown to be NP-hard, and two approximation algorithms with performance 
guarantees have been proposed. 

Interesting aspects requiring further investigation include a throughout study of the latency optimization 
problem for cyclic roadmaps, the development of more efficient approximation algorithms, and the intro- 
duction of a more general communication framework, in which the robots are allowed to communicate 
while traveling the edges of the roadmap. The study of average performance criteria and the extension 
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to dynamically changing environments are also of interest. Finally, an hardware implementation of our 
algorithms would further strengthen our theoretical findings. 
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